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Abstract — The direct-sequence code-division multiple access 
(DS-CDMA) cellular downlink is modeled by a constrained 
random spatial model involving a fixed number of base stations 
placed over a finite area with a minimum separation. The analysis 
is driven by a new closed-form expression for the conditional 
outage probability at each mobile, where the conditioning is with 
respect to the network realization. The analysis features a flexible 
channel model, accounting for path loss, Nakagami fading, and 
shadowing. By generating many random networks and applying 
a given resource allocation policy, the distribution of the rates 
provided to each user is obtained. In addition to determining 
the average rate, the analysis can determine the transmission 
capacity of the network and can characterize fairness in terms 
of the fraction of users that achieve a specified rate. The analysis 
is used to compare a rate-control policy against a power-control 
policy and investigate the influence of the minimum base-station 
separation. 

I. Introduction 

The classical approach to analyzing the cellular downlink 
is to model the network as a lattice or regular grid of 
base stations [1|. By using the geometry of the grid, which 
typically assumes a honeycomb-like arrangement of hexagons, 
along with models for the underlying propagation conditions 
and channel reuse strategies, performance metrics can be 
computed at various potential mobile locations. Often, the 
analysis focuses on the worst-case locations, which are at the 
cell edges. The hexagonal grid model was used to analyze 
the other-cell interference (OCI) of a power-controlled direct- 
sequence code-division multiple access (DS-CDMA) downlink 
in J2[ . Although conceptually simple and locally tractable, 
the grid assumption is a poor model for actual base-station 
deployments, which cannot assume a regular grid structure 
due to a variety of regulatory and physical constraints. 

Recently, a new approach to the analysis of cellular net- 
works has been proposed in [3| and [4|, which models the 
base-station locations as a realization of a random point 
process, thereby allowing the use of analytical tools from 
stochastic geometry [5|. This approach can be used to de- 
termine the performance of a typical mobile user. By com- 
bining the effects of random base-station location, fading, and 
shadowing into a single random variable, performance metrics 
such as coverage probability (the complement of the outage 

M.C. Valenti and S. Talarico were sponsored by the National Science 
Foundation under Award No. CNS-0750821 and by the United States Army 
Research Laboratory under Contract W91 1NF-10-0109. 



probability) and average achievable rate can be determined. 
Under certain limitations, the performance metrics can be 
found in closed form with surprisingly simple expressions, 
and these expressions can provide insight into the influence 
of key network and channel parameters such as the path- 
loss exponent, the density of base stations, and the minimum 
required signal-to-interference-and-noise ratio (SINR) required 
to achieve acceptable coverage. More recently, this approach 
has been extended to account for multi-tier heterogenous 
networks J6). 

The grid-based and random-spatial approaches represent 
two extreme perspectives in modeling cellular networks. 
Whereas the grid-based approach overly constrains the base- 
station locations, a pure random-spatial approach is similarly 
unrealistic because it does not place enough constraints on 
the base-station placements. For instance, the spatial model 
in Q assumes that the base stations are drawn from a two- 
dimensional Poisson point process (PPP) and that the network 
extends infinitely on the Euclidian plane. This is a poor model 
for several reasons. The first is that no network has an infinite 
area. The second is that any realization of a PPP could 
have a large, even infinite, number of base stations placed 
within a finite area. Finally, the pure PPP model does not 
permit a minimum separation between base stations, which is 
characteristic of actual macro-cellular deployments. 

In a typical modern cellular system, each base station is 
allocated a maximum total transmit power, which must be 
shared among the users in its cell through some resource 
allocation policy. Simple spatial models do not adequately 
capture the nuances of power allocation policies. For instance, 
even a simple equal-share power allocation policy is difficult to 
model analytically because the power allocated to a given user 
will depend on the number of mobiles in the same cell. As a 
complement to power control, systems often use rate control to 
ensure that the rate provided to each user is maximized under 
a constraint on outage probability. While average achievable 
rate is considered in |3], it is computed under the assumption 
that the rate provided to each mobile perfectly adapts to the 
instantaneous SINR such that the outage is zero. This overly 
optimistic assumption is not realistic for current rate-control 
implementations, which adapt the modulation and coding 
scheme to provide a particular outage probability (typically 
0.1) when averaged over the fading, but conditioned on the 
network realization. 



In this paper, we use a constrained random spatial approach 
to model the DS-CDMA downlink. The spatial model places a 
fixed number of base stations within a region of finite extent. 
The model enforces a minimum separation among the base 
stations. The model for base-station placement is a binomial 
point process (BPP) with repulsion, which we call a uniform 
clustering model. To facilitate the study of resource allocation 
policies that depend on the mobile locations, the mobiles are 
also placed according to a uniform clustering process with a 
higher density and smaller minimum separation than that used 
to place the base stations. 

The analysis in this paper is driven by a new closed- 
form expression, recently published in |7|, for the conditional 
outage probability at each mobile, where the conditioning is 
with respect to the network realization. The approach involves 
drawing realizations of the network according to the desired 
spatial and shadowing models, and then computing the outage 
probability at each realized mobile location. Because the 
outage probability at each mobile is averaged over the fading, 
it can be found in closed form with no need to simulate 
the corresponding channels. A Nakagami-m fading model is 
assumed, which models a wide class of channels, and the 
Nakagami-m fading parameters do not need to be identical 
for all communication links. This is a useful feature that can 
be used to model situations where the base station serving a 
mobile is in the line of sight (LOS) while the interfering base 
stations are non-LOS. 

By averaging over many network realizations, the mean 
outage probability can be found for a typical mobile. However, 
characterizing only the average performance is of limited 
utility, and the approach presented in this paper allows the 
analysis to extend beyond the ergodic performance of the 
typical mobile user. For instance, by averaging over the mo- 
biles that are farthest from their base station, the performance 
of cell-edge mobiles can be determined. More generally, the 
outage probability of each mobile can be constrained, and the 
statistics of the rate provided to each user can be determined 
under various power-control and rate-control policies. By plot- 
ting the complementary cumulative distribution function (ccdf) 
of the per-user rate, the variability in mobile performance can 
be visualized. 

The remainder of this paper is organized as follows. Sec- 
tion [II] presents a model of the network culminating in an 
expression for SINR. Section [In] provides an expression for 
the outage probability based on the analysis published in 
0. Section IV discusses policies for rate control and power 



control. A performance analysis is given in Section [Vj which 
compares power control with rate control on the basis of 
average rate, transmission capacity, and fairness. The Section 
also investigates the influence of the minimum base-station 



separation. Finally, the paper concludes in Section VI 
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II. Network Model 

network comprises M cellular base stations 
Xm} and K mobiles {Yi,...,Yk} placed on a 




Fig. 1. Close-up of an example network topology. Base stations are 
represented by large rilled circles, and mobiles by small dots. Cell boundaries 
are indicated, and the minimum base-station separation is rj, s = 0.25. The 
average cell load is K/M = 16 mobiles. 

Xi represents both the i fh base station and its location, and 
similarly, Yj represents the j th mobile and its location. An 
exclusion zone of radius surrounds each base station, 
and no other base stations are allowed within this zone. 
Similarly, an exclusion zone of radius r m surrounds each 
mobile, and no other mobiles are allowed within a placed 
mobile's exclusion zone. Fig.[T]shows a portion of an example 
network with average number of mobiles per cell K/M = 16, 
a base-station exclusion radius = 0.25, and a mobile 
exclusion radius r m = 0.01. The base station locations are 
given by the large filled circles, and the mobiles are dots. A 
Voronoi tessellation shows the cell boundaries that occur in 
the absence of shadowing. 

Each mobile connects to at most one base station. Let 3^ be 
the set of mobiles connected to base station Xi, and K. L — |3^| 
be the number of mobiles served by Xi. Furthermore, let g(j) 
be a function that returns the index of the base station serving 
Yj so that Yj 6 3^ if g(j) = i- If Yj cannot connect to any 
base station, which is possible when a cell runs out of available 
channels, then g(j) = 0. 

The downlink signals use orthogonal DS-CDMA sequences 
with common spreading factor G. Because the sequences 
transmitted by a particular base station are orthogonal, the 
only source of intracell interference is due to multipath and 
the corresponding loss of orthogonality. However, if the ratio 
of the maximum-power and minimum-power multipath com- 
ponents is sufficiently small (e.g., less than about 0.1G), then 
the multipath components will have negligible effect. For this 
reason, we neglect the intracell interference and assume that 
intercell interference is the only source of interference. 

The signal is transmitted by base station Xi to mobile Yj 



and area A r 



The variable with average power Pi _j. We assume that the base stations 



transmit with a common power P such that 
1 
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for each i, where f p is the fraction of the base-station power 
reserved for pilot signals needed for synchronization and 
channel estimation. Power allocation strategies are considered 
later in this paper. 

Using spreading sequences with a spreading factor of G 
directly reduces the power of the intercell interference. While 
the intracell sequences transmitted by a particular cell's base 
station are synchronous, the varying propagation delays from 
the other base stations cause the intercell interference to be 
asynchronous. Because of this asynchronism, the intercell 
interference is further reduced by the chip factor /i(rj j.), which 
is a function of the chip waveform and the timing offset n t k at 
the mobile between the signal received from interfering base 
station Xk and the signal received from serving base station 
Xi. When the timing offset Ti,u is assumed to have a uniform 
distribution over the chip interval, then the expected value of 
h(i~i,k) is 2/3 (8). ^ is assumed henceforth that G/h(ji t k) is 
a constant equal to G/h at each mobile in the network. 

The power of Xi received at Yj also depends on the 
fading and path loss models. We assume that path loss has 
a power-law dependence on distance and is perturbed by 
shadowing. When accounting for fading and path loss, the 
despread instantaneous power of X{ at mobile Yj is 
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where gi j is the power gain due to fading, £jj is a shad- 
owing factor, and /(•) is a path-loss function. The {gij} 
are independent with unit-mean, and gij = afj, where 
Ojj is Nakagami with parameter rriij. While the {gij} are 
independent from mobile to mobile, they are not necessarily 
identically distributed, and each mobile can have a distinct 
Nakagami parameter rriij. When the channel between X; and 
Yj experiences Rayeigh fading, rrii j — 1 and <j>.y is exponen- 
tially distributed. In the presence of log-normal shadowing, 
the are i.i.d. zero-mean Gaussian with variance a].. In 

the absence of shadowing, = 0. For d > do, the path-loss 
function is expressed as the attenuation power law 



f(d) = 



(3) 



where a > 2 is the attenuation power-law exponent, and that 
do is sufficiently large that the signals are in the far field. 

The base station X g (j) that serves mobile Yj is selected to 
be the one with index 

g(j) = argmax 



{lO^/MfiWXt-YjW)} 



(4) 



which is the base station with minimum path loss to Yj. In the 
absence of shadowing, it will be the base station that is closest 



to Yj, In the presence of shadowing, a mobile may actually 
be associated with a base station that is more distant than the 
closest one, if the shadowing conditions are sufficiently better. 
We assume a maximum of G orthogonal spreading sequences 
per cell, and once G users are connected to a base station, no 
more can be served. 



The instantaneous SINR at mobile K, is 
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where JV is the noise power. Substituting Q and Q into |5]) 
yields 
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where T = dyPo/J\f is the signal-to-noise ratio (SNR) at a 
mobile located at unit distance when fading and shadowing 
are absent, and 

fiu = ^io«-/ 10 ||x-r,||- Q (7) 

is the normalized power of Xi at receiver Yj before despread- 
ing. 

III. Outage Probability 
Let f3j denote the minimum SINR required by Yj for 



reliable reception and flj 



1,3) • 



, Omj} represent the 



set of normalized despread base-station powers received by 
Yj. An outage occurs when the SINR falls below f)j. As 
discussed subsequently, there is a relationship between the 
SINR threshold and the supported rate of the transmission. 
Conditioning on Slj, the outage probability of mobile Yj is 



ej = P[ 7j <l3j\nj]. 



(8) 



Because it is conditioned on fij, the outage probability 
depends on the particular network realization, which has 
dynamics over timescales that are much slower than the fading. 
By defining a variable 
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the conditional outage probability may be expressed as 

ej = P[Z 3 <Y- 1 \n j \=F Z] {T- l \Slj) 



(10) 



which is the cumulative distribution function (cdf) of Zj 
conditioned on ttj and evaluated at r _1 . 

Define Fz^z^j) = 1 — Fz^zl^lj) to be the complemen- 
tary cdf of Zj conditioned on S~2j . Restricting the Nakagami 
parameter rrig^ j between mobile Yj and its serving base 



station X g (j\ to be integer-valued, the complementary cdf of 
Zj conditioned on Slj is proved in Q to be 
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where too — m g(j).j> A) = P^o/^o, 
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the summation in ( fl~3] > is over all sets of indices that sum to 
fc, and 
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IV. Policies 



rriii+e 



(14) 



A key consideration in the operation of the network is the 
manner that the total power Pq transmitted by a base station 
is shared by the mobiles it serves, which influences the rates 
provided to each user. This section discusses two options for 
allocating rate and power, rate control and power control. 

A. Rate Control 

A simple and efficient way to allocate Pq is with an 
equal-share policy, which involves base station Xi,i = g(j), 
transmitting to mobile Yj with power 

Po 
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Under this policy, the SINR will vary dramatically at each 
mobile. If a common SINR threshold (3j is used by all mobiles, 
then the outage probability will likewise be highly variable. 
Instead of using a fixed threshold, the threshold ftj of mobile 
Yj can be selected such that the outage probability of mobile 
Yj satisfies the constraint tj = e. A constraint of i = 0.1 is 
typical and appropriate for modern systems that use a hybrid 
automatic repeat request (HARQ) protocol. 

For a given j3j, there is a corresponding transmission rate 
Rj that can be supported. Let Rj = C (j3j ) represent the rela- 
tionship between Rj, expressed in units of bits per channel use 
(bpcu), and j3j. The relationship depends on the modulation 
and coding schemes used, and typically only a discrete set of 
Rj can be supported. While the exact dependence of Rj on 
f3j can be determined empirically through tests or simulation, 
we make the simplifying assumption when computing our 
numerical results that C(j3j) = log 2 (l + fij) corresponding 
to the Shannon capacity. This assumption is fairly accurate 
for capacity-approaching codes with an infinite number of 
possible rates. This is a reasonable model for modern cellular 
systems, which use turbo codes with a large number of 
available rates. 



With an equal power share, the number of mobiles JQ in 
the cell is first determined, and then the power share given 
to each is found from (JT3J. For each mobile in the cell, the 
corresponding (3j that achieves the outage constraint tj — t is 
found by inverting the outage probability expression given in 



Section III Once f3j is found, the corresponding Rj is found 
by using the function Rj = C(j3j). The rate Rj is adapted by 
changing the number of information bits per channel symbol. 
The processing gain G and symbol rate are held constant, so 
there is no change in bandwidth. Because this policy involves 
fixing the transmit power and then determining the rate for 
each mobile that satisfies the outage constraint, we refer to it 
as rate control. 

B. Power Control 

A major drawback with rate control is that the rates provided 
to the different users in the network will vary significantly. 
This variation will result in unfairness to some users, particu- 
larly those located at the edges of the cells. To ensure fairness, 
Rj could be constrained to be the same for all Xj £ 3^. The 
power transmitted to mobile Yj is then found such that the 
mobile has outage probability tj = e, while the total power 
provided to all mobiles in the cell satisfies power constraint 
([TJ. Because the power transmitted to each mobile is varied 
while holding the rate constant for all mobiles in the cell, we 
refer to this policy as power control. Note that, while the rate 
is the same for all users within a given cell, it may vary from 
cell to cell. In particular, the rate of a given cell is found by 
determining the value of R that allows the outage constraint 
tj = t and power constraint ([T]i to be simultaneously met. 

V. Performance Analysis 

Under outage constraint e, the performance of a given 
network realization is largely determined by the set of achieved 
rates {Rj} of the K users in the network. Because the network 
realization is random, it follows that the set of rates is also 
random. Let the random variable R represent the rate of 
an arbitrary user. The statistics of R can be found for a 
given class of networks using a Monte Carlo approach as 
follows. Draw a realization of the network by placing M 
base stations and K mobiles within the disk of radius r net 
according to the uniform clustering model with minimum 
base-station separation rf, s and minimum mobile separation 
r m . Compute the path loss from each base station to each 
mobile, applying randomly generated shadowing factors if 
shadowing is present. Determine the set of mobiles associated 
with each base station. At each base station, apply the power 
allocation policy to determine the power it transmits to each 
mobile that it serves. By setting the outage equal to the 
outage constraint, invert the outage probability expression to 
determine the SINR threshold for each mobile in the cell. By 
applying the function Rj = C(jj), find the rate of the mobile. 
Repeat this process for a large number of networks, all with 
the same spatial constraints. 

Let E[R] represent the mean value of the variable R, 
which can be found by numerically averaging the values of 
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Fig. 2. Transmission capacity of the Example as a function of K/M with 
rate control and power control. 

R obtained using the procedure described in the previous 
paragraph. While E[R] is a useful metric, it does not account 
for the loss in throughput due to the inability to successfully 
decode during an outage, and it does not account for the spatial 
density of transmissions. These considerations are taken into 
account by the transmission capacity, defined as J9) 

r = X(l-e)E[R] (16) 

where A = K/A net is the density of transmissions in the 
network. Transmission capacity can be interpreted as the 
spatial intensity of transmissions; i.e., the rate of successful 
data transmission per unit area. 

As an example, consider a network with M = 50 base 
stations placed in a network of radius r net = 2 with base- 
station exclusion zones of radius = 0.25. A variable 
number K of mobiles are placed within the network using 
exclusion zones of radius r m — 0.01. The outage constraint 
is set to e = 0.1 and both power control and rate control are 
considered. A mobile in an overloaded cell is denied service. 
In particular, the Ki — G mobiles whose path losses from the 
base station are greatest are denied service, in which case they 
do not appear in the set J^i for any i, and their rates are set 
to Rj = 0. The SNR is set to F = 10 dB, the fraction of 
power devoted to pilots is f p = 0.1, and the spreading factor 
set to G = 16 with chip factor h = 2/3. The propagation 
environment is characterized by a path-loss exponent a = 3, 
and the Nakagami factors are rrii.j = 3 for i — g(j) while 
rriij = 1 for i ^ g(j); i.e., the signal from the serving base 
station experiences milder fading than the signals from the 
interfering base stations. This is a realistic model because the 
signal from the serving cell is likely to be in the line of sight 
(LOS), while the interfering cells are typically not LOS. 

Fig. [2] shows the transmission capacity, as a function of the 
ratio K/M, of rate control and power control in an unshad- 
owed environment as well as in the presence of shadowing 
with a s = 8 dB. The figure shows that the transmission 
capacity under rate control is higher than it is under power 
control. This disparity occurs because mobiles that are close 
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Fig. 3. Average rate of the Example as a function of K/M in the presence of 
shadowing with rate control and power control. For rate control, the averaging 
is done over all mobiles and over just the cell-edge mobiles. With power 
control, all mobiles in a cell are given the same rate. 

to the base station are able to be allocated extremely high rates 
under rate control, while with power control, mobiles close to 
the base station must be allocated the same rate as mobiles at 
the edge of the cell. As the network becomes denser (K/M 
increases), shadowing actually improves the performance with 
rate control, while it degrades the performance with power 
control. This is because shadowing can sometimes cause the 
signal power of the base station serving a mobile to increase, 
while the powers of the interferers are reduced. The effect 
of favorable shadowing is equivalent to the mobile being 
located closer to its serving base station. When this occurs with 
rate control, the rate is increased, sometimes by a very large 
amount. While shadowing does not induce extremely favorable 
conditions very often, when it does the improvement in rate 
is significant enough to cause the average to increase. On the 
other hand, a single user with favorable shadowing conditions 
operating under power control will not improve the average 
rate because all mobile in a cell receive the same rate. 

While rate control offers a higher average rate than power 
control, the rates it offers are much more variable. This 
behavior can be seen in Fig. [3] which compares the rates of 
all users against those located at the cell edges. In particular, 
the figure shows the rate averaged across all mobiles for both 
power control and rate control, as well as the rate averaged 
across just the cell-edge mobiles for rate control, where the 
cell-edge mobiles are defined to be the 5 percent of the mobiles 
that are furthest from their serving base station. The average 
rate of cell-edge mobiles is not shown for power control 
because each cell-edge mobile has the same rate as all the 
mobiles in the same cell. As seen in the figure, the performance 
of cell-edge mobiles is worse with rate control than it is with 
power control. 

The fairness of a particular power allocation policy can 
be further visualized by plotting the ccdf of R, which is the 
probability that R exceeds a threshold r; i.e. P[R > r]. Fig. [4] 
shows the ccdf of R for the Example with shadowing (a s = 8 




Fig. 4. Ccdf of R with either rate control or power control for a lightly 
loaded system (K/M = 4) and a moderately-loaded system (K/M = 12). 




Fig. 5. Transmission capacity for rate and power control as function of rf, a 
for K/M = 8 in Mixed Fading and Shadowing (cr s = 8 dB). 

dB) and either rate control or power control. Two system 
loads are considered: a lightly loaded system (K/M = 4) 
and a moderately-loaded system (K/M = 12). The ccdf 
curves for power control are steeper than they are for rate 
control, indicating less variability in the provided rates. The 
lower variability in rate corresponds to improved fairness. For 
instant, with a load K/M = 4, almost all (99.9%) users are 
provided with rates of at least r = 0.5 under power control, 
while with rate control 96% of the users are provided with 
rates of at least r — 0.5 implying that a significant fraction of 
the users are provided with lower rates. 

A key feature of the proposed analysis is that it permits 
a minimum spacing of r^s around each base station. The 
dependence of the transmission capacity on is shown in 
Fig. [5] for both rate control and power control with a system 
load of K/M — 8. Curves are shown for two values of path- 



loss exponent: a = 3 and a = 4. Except for the values of 
and a, the parameters are the same as in the Example with 
shadowing. Increasing rj, s improves the transmission capacity. 
The improvement is slightly more pronounced for rate control 
than for power control, and is more pronounced for the higher 
value of a. 

VI. Conclusion 

This paper has presented a powerful new approach for 
modeling and analyzing the cellular downlink. Unlike simple 
spatial models, the model in this paper allows constraints to be 
placed on the distance between base stations, the geographic 
footprint of the network, and the number of base stations 
and mobiles. The analysis features a flexible channel model, 
accounting for path loss, shadowing, and Nakagami-m fading 
with non-identical parameters. The proposed analytical ap- 
proach provides a way to compare various access and resource 
allocation techniques. As a specific application of the model, 
the paper models a direct-sequence CDMA network, analyzes 
it using realistic network parameters, and compares the per- 
formance of two resource allocation policies. While this paper 
only considers resource allocation policies that implement 
either rate control or power control, other policies that fall 
between these two extremes could be considered. While in 
this paper mobiles in overloaded cells were denied service, 
the work can be extended to analyze reselection schemes, 
whereby mobiles in overloaded cells attempt to connect to 
another nearby base station serving an underloaded cell. This 
work could be extended to analyze the uplink and to model 
other types of access, such as orthogonal frequency-division 
multiple access (OFDMA). It could furthermore be extended 
to handle sectorized cells and coordinated multipoint strategies 
involving transmissions from multiple base stations. 
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